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We investigate the distribution of number of photons emitted by a single molecule undergoing a 
spectral diffusion process and interacting with a continuous wave field. Using a generating function 
formalism an exact analytical formula for Mandel's Q parameter is obtained. The solution exhibits 
transitions between; (i) Quantum sub-Poissonian and Classical super-Poissonian behaviors, and (ii) 
fast to slow modulation limits. Our solution yields the conditions on the magnitude of the spectral 
diffusion time scales on which these transitions are observed. We show how to choose the Rabi 
frequency in such a way that the Quantum sub-Poissonian nature of the emission process becomes 
strongest, we find ^^ = where Fsd (P) is the spectral diffusion (radiative) contribution to 

the width of the line shape. 



PACS numbers: 82.37.-j, 05.10.Gg, 33.80.-b, 42.50.Ar 
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Physical, Chemical, and Biological systems are inves- 
tigated in many laboratories using single molecule spec- 
troscopy U . The investigation of the distribution of num- 
ber of photons emitted from a single molecule source is 
the topic of extensive theoretical research e.g. [1,11,^13' 
An important mechanism responsible for fluctuations in 
photon statistics from single molecule sources is spectral 
diffusion. In many cases the absorption frequency of the 
molecule will randomly change due to different types of 
interactions between the molecule and its environment 
(e.g- [ElEQilS and Ref. therein). For example in low 
temperature glasses flipping two level systems embedded 
in the glassy environment, induce spectral jumps in the 
absorption frequency of the single molecule under inves- 
tigation In this way the molecule may come in and 
out of resonance with the continuous wave laser field with 
which it is interacting. 

Obviously a second mechanism responsible for fluctu- 
ations of photon counts is the quantum behavior of the 
spontaneous emission process [2, Il0| | . In his fundamen- 
tal work Mandel showed that a single atom in the 
process of resonance fluorescence, in the absence of spec- 
tral diffusion, exhibits sub-Poissonian photon statistics 
[T^ . Photon statistics is characterized by Mandel's Q 
parameter 



N 



- 1 



(1) 



where N is the number of emitted photons within a 
certain time interval. The case Q < is called sub- 
Poissonian behavior, while Q > is called super- 
Poissonian. Sub-Poissonian statistics has no classical 
analog [13. Briefly, the effect is related to anti-bunching 
of photons emitted from a single source and to Rabi- 
oscillations of the excited state population which favors 
an emission process with some periodicity in time. Sub- 
Poissonian statistics and anti-bunching were detected in 
several single molecule experiments, for example J-S, 14] . 

In this Letter we obtain an exact analytical expression 
for the Q parameter in the long time limit, for a single 



molecule undergoing a stochastic spectral diffusion pro- 
cess. For that aim we use the powerful generating func- 
tion method of Zheng and Brown If the spectral 
diffusion process is slow enough super-Poissonian behav- 
ior of Q is expected (see details in text). Our analytical 
expressions classify the transitions between sub and su- 
per Poissonian statistics. They give the conditions on the 
spectral diffusion time scale and magnitude of spectral 
jumps for sub-Poissonian non-classical behavior to be ob- 
served. Our result is valid for weak and strong excitation 
(i.e arbitrary Rabi frequency). It yields the lower bound 
on Q. The solution shows how in experiment we may 
choose the Rabi frequency so that the quantum nature 
of the photon emission process becomes larger, namely 
how to minimize Q in the Sub-Poissonian regime. This 
is important for the efficient detection of quantum effects 
in single molecule spectroscopy, since choosing too small 
or too large values of the Rabi frequency results in very 
small and hence undetectable values of Q. Finally our 
result is used to test the semi-classical linear response 
theory (i.e^ weak Rabi frequency) of Barkai, Jung, and 
Silbey [a |8|. The semi-classical theory yields Q > 0, 
while the main focus of this manuscript is on the quan- 
tum regime Q < 0. 

Our starting point are the Zheng-Brown (l5| general- 
ized optical Bloch equations describing a chromophore 
with single excited and ground state 

U (s) = -^U (s) + 6L{t)V (s) 

V (s) = -dL{t)u (s) - §v (s) - nw (s) 



w (s) = nv (s) - ^ (1 s) w (s) - § (1 + s) y (s) 



-L^l-s)Wis)-^il-s)y{s). 



(2) 

These equations are exact within the rotating wave ap- 
proximation and optical Bloch equation formalism. In 
Eq. |(5J) r is the spontaneous emission rate of the elec- 
tronic transition and is the Rabi frequency. The time 
evolving detuning is SL(t) = ujl — ujq — Auj{t), where 
LUL (ij^o) is the laser frequency (the molecule's bare fre- 
quency), and Auj{t) is the stochastic spectral diffusion 
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process. The information about the photon statistics is 
contained in the moment generating function (15| 

N=0 

where PN{t) is the probabiUty of N emission events in 
the time interval (0,i). Thus the solution of Eq. Q 
yields the moments of N 

]v(I)-2y(i) Wit) ^2y"{i) + 2y{i), (4) 

with which the Mandel Q parameter can in principle be 
obtained. In Eq. (0J, and in what follows, we use the no- 
tation ■^g{s)\s=i = and similarly for second order 
derivatives with respect to s. The physical meaning of 
U{s), V(s), and yV(s) and their relation to the standard 
Bloch equation was given in and in closely related 
work by Mukamel [T^ , some discussion on this issue will 
follow Eq. © . Note that when s — > 1 the damping terms 
in Eq. become small [i.e. the (1 — s)r/2 terms], hence 
relaxation of the generalized Bloch equations in the im- 
portant limit of s — > 1 is slow. 

In what follows we will consider the moments 
N{t), N-^it). For this aim it is useful to de- 
rive equations of motion for the vector z — 

{w(i),v(i),w(i),3^(i),w'(i),v'(i),w'(i),y(i),y'(i)}. 

Taking the first and the second derivative of Eq. 
with respect to s and setting s = 1, we find 

z^M{t)z (5) 

where M{t) is a 9 x 9 matrix M{t) = 
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The first three lines of M(t) describe the evolution of 
Z//(l),V(l),yV(l), these are the standard optical Bloch 
equations in the rotating wave approximation. The 
fourth line of M(i) is zero, it yields 3^(1) = 0, this equa- 
tion describes the normalization condition of the problem 
namely 3^(1) = 1/2 for all times t [to see this use Eq. (O 
and X]w=o^^(^) ~ The evolution of the remaining 
terms Z^'(l), V'(l), W(l), y (1), y (1) are of current in- 
terest since they describe the fluctuation of the photon 
emission process. Solutions of time dependent equations 
like Eq. Q are generally extremely difficult to obtain, 
a formal solution is given in terms of the time ordering 
operator T, = rexp[/g M(t)dt]z(0). 

Eq. (O yields a general method for the calculation of 
Q for a single molecule undergoing a spectral diffusion 
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FIG. 1: In (a) we show the line shape litoi,) vs ujl for F = 40 
MHz, Q. = r/V2, and v = 2V. The line exhibits the well 
known transition between fast (_R = lOOF dotted solid curve) 
and slow modulation limits (-R = F/IO solid curve). In (b) 
we show the Q parameter, using the same parameters. Three 
type of behaviors are observed: (i) for fast modulation the Q 
parameter attains its minimum on zero detuning, (ii) For an 
intermediate modulation, ii = r/2 (dotted curve) we obtain 
Q < however now the minimum of Q is on li^l — li^o = 
±v. (iii) For the slow process Q exhibits super-Poissonian 
statistics Q > 0. Note that we used Q, = r/\/2 in order to 
maximize the absolute value of Q in the fast modulation limit. 



process. The aim of this Letter is to obtain an exact so- 
lution of the problem for an important stochastic process 
used by Kubo and Anderson to investigate characteris- 
tic behaviors of line shapes. We assume Aw(f) = vhit) 
where v describe frequency shifts, and hit) describes a 
random telegraph process: hit) — 1 or hit) = —1. The 
transition rate between state up and state down (-) 
and vice versa is R. This stochastic process was used to 
describe behavior of line shapes in many systems ;f|, , 
here our aim is to calculate Q describing the line shape 
fluctuations. 

We use Burshtein's method 0, of marginal aver- 
ages, to solve the stochastic differential matrix equation 
(O. The method yields the average {z) with respect to 
the stochastic process. We will calculate {z) in the limit 
of long times, and then obtain the steady state behavior 
of the line shape / (wl) = limt^oo jt i-^it)) and the Man- 
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del Q parameter. Let {z)± be the average of z{t) under 
the condition that at time t the value oih{t) = ±1 respec- 
tively. {z)± are called marginal averages, the complete 
average is (z) — (2;)+ + (z)-. The equation of motion for 
the marginal averages is 



M+ - RI 
RI 



RI 

- RI 



(z) 



(7) 



In Eq. Q the matrix M± is identical to matrix M in 
Eq. (jSJ when replaced by = ujl — ^^o T t^j 

and / is a 9 X 9 identity matrix. The 18 x 18 ma- 
trix in Eq. O is traceless, hence the time depen- 
dent solution of Eq. Q is not obtained by direct in- 
version of it. We found the long time behavior of 
Eq. iQ after a long journey which involved the fol- 
lowing five steps. 1. We obtained the steady state 



solution 
2. We 

step 



where 



A = 



for ({ZY(l)±,V(l)±,W(l)±,3^(l)±,y(l)±}). 
used the solutions obtained in previous 
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and b{t) = (0, 0, &+(<), 0, 0, with 



(9) 



b±{t) = - ^W^'{1)± - X + ^ + ± ^ [W^^(l)- - W^^(l)+] . (10) 

and 

- I [(^"')33 + (^") J W^^(l)- = I [(^")63 + (^")66] • (11) 

We note that W''{1)± yields the steady state (ss) marginal averages of the population difference between the upper 
excited state and ground state. 3. Using Mathematica we found an analytical expression for A^^, the inverse of the 
Matrix A, in terms of the parameters of the problem R, ly, T, lol — ^0, 4. Using A~^ and Eq. (jSJ we find the long 
time behavior of the vector x. We then use Eq. Q to obtain steady state behavior of (3^"(1)*''')±. 5. Finally we 
obtain the Q parameter using: Q = [(y'(l)"") - 2{y'(lY^)'^] /{y'{lY^). 
Following these steps we obtain the main result of this Letter 

Q ^ Y E E (^"'^"')., + jnrr\ [w''(i)+ - w''(i)-]' + 77^ E [>^''(i)'= (1 + 4W^'(i)'=)] - (12) 

1=3,6 j=3, 6 ^ ^ ^' k=± 



where the line is 



(13) 



The general expression for Q in terms of original param- 
eters of the problem can be found in Below we 
will analyze the physical behaviors of Q obtaining simple 
equations in limiting cases. 

In Fig. n we show the line shape and the Q param- 
eter. We choose = r/\/2, fix and vary the rate 
R. The line exhibits two well known behaviors: a fast 
modulation limit where the line is Lorentzian and a slow 
modulation limit where the line exhibits two peaks on 
ldl — ±J^. The Q parameter exhibits three behaviors: 
(i) fast modulation limit, where the single molecule ex- 
hibits a sub-Poissonian behavior, and Q attains its min- 
imum on zero detuning, (ii) Intermediate modulation 
limit where photon statistics is still sub-Poissonian, how- 



ever now Q attains its minimum on ljl = '^v- And (iii) a 
slow modulation limit where a super-Poissonian behavior 
is observed. We now analyze these limits in some detail. 
Without loss of generality we use units where F = 1, and 
set ojQ = hence lol is the detuning frequency. 

For zero detuning, ljl = 0, we obtain from Eq. (|12|l 



Q 



4 (-1 + 4i?) -t-3 (1 -H4i?)^ 



[Av^ + {\ + 29?) (l + 4i?)]'' 



(14) 



Eq. (|14|l indicates a transition from sub-Poissonian 
statistics {Q < 0) to super-Poissonian statistics {Q > 0) 
when 4z/2 = 3(1 + 4i?)2/(l - 4i?). For fast processes 
satisfying R > 1/4, (i.e. R > r/4 in natural units) 
we find sub-Poissonian behavior for all values of v and 
n. As mentioned in introduction an important issue is 
how to choose D, in such a way to minimize Q in the 
sub-Poissonian regime (i.e., how to maximize the quan- 
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turn nature of the photon statistics). Using Eq. H14() 
it is easy to find the lower bound Q > —3/4. The ab- 
solute minimum Q = -3/4 is found when n = 1/^2 
and V = 0, namely for a stable molecule whose ab- 
sorption frequency is not fluctuating. For fixed v and 
R the largest quantum fluctuations are observed when 
nl^.^^ = (1 + 4i/2 + 4i?)/[2(l + 4i?)]. Note however that 
for fast processes v « R, the parameters v and R can- 
not be directly deduced in experiment, due to motional 
narrowing effect. In Eq. (|17|l we will give an expression 
for rJmin in terms of physical observables. 

Consider the fast modulation limit. If we take i? — > oo 
keeping and v fixed we obtain from Eq. \i'2.l 



lim Q 

i?-»oo 



29? (3 - AujI) 



(1 + 2^2 



(15) 



This result was derived by Mandel for an atomic transi- 
tion in the absence of spectral transition. As expected 
when the time scale for spectral jumps becomes very 
short, the molecule does not respond to the stochastic 
fiuctuation. A more physically interesting case is to let 
i? — > cxD and — > oo but keep the spectral diffusion (SD) 
contribution to the line width Fgo = / R remain finite. 
In this limit, the line shape is Lorentzian and it exhibits 
motional narrowing (i.e. as R is increased the line be- 
comes narrower). In this fast limit we obtain from Eq. 



Q 



fast 



[3 - 



^ SD 



4ujL^+TsD (7 + 4c^i)] 



[l + TsD 2 + 2^(2 + 2rsD (1 + 1^2) 



(16) 

Note that photon statistics are sub-Poissonian provided 
that the detuning is not too large. The Q parameter 
obtains its minimum when 



(l+T; 



SD) 



2(i + rsD) 



(17) 



and for zero detuning Q^in = -(Esd + 3)/4(l -|- Tsd)- 
The simple Eq. H17(l is important, since it shows how to 
choose the Rabi frequency in order to obtain strong sub- 
Poissonian behavior. Since Psd is a physical observable 
(unlike R and i') Eq. (|17|l is of practical value. 



To investigate the intermediate and slow modulation 
limit we consider the value of Q for = z/. We consider 
a strong coupling limit z^ >> 1, (e.g. spectral jumps of 
the order of v ~ GHz >> F). From Eq. ((T^ we obtain 

lim Qujl=v = 



(1 + 2R) n"^ [-1 + 16i?2 + 8i?3 - 2r22 + 2i? (2 + Vi^)] 
2R[l + AR^ + 29P- + 2R{2 + 9?)f 

(18) 

It is easy to see that Eq. (I18|l exhibits both sub- 
Poissonian and super-Poissonian behaviors. When the 
process is very slow, namely i? ^ 0, we obtain 



lim Q^^=u 



fi2 



^^oo 2R{1 + 2Q?) 



(19) 



a super-Poissonian behavior. In the intermediate modu- 
lation limit, for example when R = 1, we obtain 



lim 



8ir22 



2 (9 + 4r22 



(20) 



a sub-Poissonian behavior. In this case the fluctuations 
are strongest when ri^i,, = 9/4 and then Qmin = —9/32. 
It is important to notice that in the same limit — > oo 
we have for zero detuning \Q\ oc v^^ hence in the 
slow and intermediate modulation limits the fluctuations 
on ljl — ±1^ are much stronger than the fluctuations on 
zero detuning (an example is in Fig. 

The detailed interaction between a chromophore and 
a bath is usually extremely complex, still for many fast 
processes line shapes are simple Lorentzians (universal- 
ity). We expect that certain behaviors of Q obtained in 
this work are also universal, mainly the behavior in the 
fast modulation limit and the transition between fast, in- 
termediate, and slow limits. We also expect, but have no 
proof, that the maximum of the quantum fluctuations, is 
attained close to ^^„i„/F = (Fsd +r)/2 [i.e., Eq. ^ in 
natural units and for zero detuning] , also in other models 
of line broadening. Finally, the optimization of external 
fields to obtain strong sub-Poissonian and anti-bunching 
signals from single molecules is a new topic of research 
worthy of further investigation. 
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I 

I. EXACT SOLUTION 

The exact solution for the Q parameter, with F = 1, in the long time limit, Q — Numerator[(5]/Denominator[Q] is 

Denominator [Q] = 

R {Aljl^ (l + 2i?) + 4;/2 (i + 6i? + 8i?2) + (l + 4i?) (l + 6 i? + 8 i?^ + 2 f^^)) (i6jy4 (i + 2i?) + 

8z^2 il^AijjL^ (1 + 2i?) + 21^2 + 2i? (3 + 4i?) (l + f^^)) + (l + 4cjl^ + 2 r?^) (21) 

(4wL^ (l + 2i?) + (l + 4i?) (l + 6i? + 8i?2 + 2r22)))^ 



Numerator[(5] = 

-2^22 [2hQu^ R{l + 2Rf (-l + 16i?2) - {-i + Aujl) R {lul^ (4 + 8 i?) + (1 + 4 i?) {l + 6 R + 8 R^ + 2n^)f 
-128 (l + 2i?)^ (4wL^ (l + 4i? + 4i?2 + 16i?3 + 32i?'') -3i? (l + 4i?) (24i?2 + 32i?3-r22^4i? {l + fl^))) 
+32;/4 (l + 2i?) (l6 uj + 2 Rf {2 + 5R + 8R^) +Alol^ (l + 6i? + 8i?2) (-2 + 112 i?"* + 64 i?^ - 4 f^^^ 
i?2 (76-8^22) (-4 + f72) + 16i?3 (10 + 172)) +3i? (l + 4i?) (l + 832 i?* + 512 i?^ + _ 3 f]* + 32 i?^ (l6 + 5f72) + 
4i?2 (37 + 321^2) +R (20 + 26172 + 81^-*))) -81^2 (q^ul^ {I + 2 Kf - R [l + A R) (l + 6 i? + 8 i?2 + 2 1^2)^ (4+152i?2 
+160i?3-172_^ 4^ (11 + 172)) - iQujL^{l + 2Rf (-2 + 40i?3 -4172 -4i?2 (_5 + 6172) ~i? (4 + 17 172)) - 4cjl2 
(l + 2i?) (7680i?6 + 2048 i?^ + 256 i?5 (37 + I72) - (l + 2 172)% 32 i?^ (169 + 36172) _^ g^a (191 + 108 172 - 8 17^) + 
4i?2 (47 + 43 172 ~ 12 174) - 2 i? (-1 + 5 172 + 14 17*)))) . 

(22) 

The line shape (F = 1) is 

/(^l) = 



http: / /www. aip .org. /pubservs /epaps .html 
[20] B. J. Berne, and R. Pecora Dynamic Light Scattering 
Dover Publication, Chapter 6. 



[172 (l + 32i?2 + 32i?3 + 4j^2 (i + 6i? + 8i?2) +2172 + 4cj2 +2i? (5 + 4172 + 4cj2))] / 
[l&v^ (l + 2i?) + 8z/2 (i + 2172 + 8i?2 (1 + 172) -4t^2 (6 + 6 172 - 8 c^i )) + 
(l + 2172 + 4tj2) (i + 32i?2_^32i?3^2172 + 4cj2 +2i? (5 + 4172 + 4cj2))] 



(23) 



